Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN KTpX style file v2.2) 



On the diversity of short GRBs 



3 

On 



Stephan Rosswog 1 and Enrico Ramirez-Ruiz 2 

1. Department of Physics and Astronomy, University of Leicester, LEI 7RH, Leicester, UK. 

2. Institute of Astronomy, Madingley Road, Cambridge, CB3 0HA, UK. 

rn ' 
O 
O 
(N 

ABSTRACT 

Hydrodynamical simulations of the last inspiral stages and the final coalescence of a double 
neutron star system are used to investigate the power of the neutrino-driven wind, the energy 
and momentum of the fireball produced via ^P-annihilation, and the intensity and character 
of their interaction. It is argued that the outflow that derives from the debris will have enough 
J> ■ pressure to collimate the relativistic fireball that it surrounds. The low luminosity relativistic 

' jet will then appear brighter to an observer within the beam although most of the energy of 

the event is in the unseen, less collimated and slower wind. This model leads to a simple 
physical interpretation of the isotropic luminosities implied for short GRBs at cosmological 
■ distances. A wide variety of burst phenomenology could be attributable to the dependence of 

£C) ' the neutrino luminosity on the initial mass of the double NS binary. 
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1 INTRODUCTION 



Gamma-ray bursts (GRBs) come in two basic types with short (< 2 s) and long (> 2 s) duration (Kouveliotou et al. 1993). So far, only the 
. latter class has been relatively well-studied - several long bursts have been associated with afterglows that have been observed from radio 
frequencies to X-ray energies (see Meszaros 2002 for a recent review). By contrast recent X-ray, optical and radio searches have not yet been 
successful for short GRBs (although see Lazzati, Ramirez-Ruiz & Ghisellini 2001). These bursts - which make up about one-third of all 
observed GRBs - differ markedly from the long ones not only in duration but also in having a larger fraction of high-energy gamma rays in 
their energy distribution. The most widely discussed possibility is that they result either from the merger of two compact neutron stars or of a 
neutron star (NS) and a black hole (BH; Lattimer & Schramm 1976; Paczynski 1986, 1991; Eichler at al. 1989; Narayan, Paczynski & Piran 
1992; Mochkovitch et al. 1993; Kluzniak & Lee 1998; Rosswog et al. 1999; Ruffert & Janka 1999; Lee & Ramirez-Ruiz 2002; Rosswog & 
Davies 2002). 

The merger of two neutron stars produces a remnant with huge energy reservoirs. The difficult question to answer is how gravity and 
spin can conspire to form the outflowing relativistic plasma that is required to explain the observations. The conventional view is that the 
released energy is quickly and continuously transformed into a radiation-dominated fluid with a high entropy per baryon (Cavallo & Rees 
1978). This fireball is then collimated into a pair of anti-parallel jets. There is as yet no universally accepted explanation for jet collimation. 
The difficulty of the problem is that there are interactions between the central object (black hole or other), the inflow (i.e. disk), the outflow 
(.i.e.wind) and the jet, and most of the dissension comes about assessing condition and strength of these interrelationships (e.g. Blandford 
2002). Particular importance is attached here to the interaction of the jet with the less collimated and slower outflow associated with the 
disk. The most prominent feature of such interaction is that the wind may be responsible for collimating the jet (Levinson & Eichler 2000; 
hereafter LE), and it is to this problem that we have turned our attention. 

We used three-dimensional high-resolution calculations of NS coalescences to study the neutrino emission from the hot merger remnant 
of a NS binary encounter. The issues we investigate here include the properties of the neutrino-driven outflow (from the debris), the energy 
and momentum injected into the region above the remnant via vv annihilation, and the intensity and character of their interaction. We estimate 
the dependence of the neutrino luminosity on the initial mass of the double NS binary and discuss the possible variety of afterglow behaviour 
that is expected from fireballs that are collimated by a surrounding baryonic wind. Our findings may be useful for designing search strategies 
for detecting the afterglows of short GRBs. 



2 Rosswog & Ramirez-Ruiz 



2 NEUTRINO EMISSION IN NS COALESCENCE 

We have recently performed a series of three dimensional, high-resolution simulations of the last inspiral stages and final coalescence of a 
double NS using the smoothed particle hydrodynamics method with up to « 10 6 particles (the reader is referred to Rosswog & Davies 2002 
for further details on the hydrodynamic evolution). Aware of its decisive role in the thermo-dynamical evolution of the debris, we have used a 
realistic equation of state for hot, dense nuclear matter, which is based on the tables provided by (Shen et al. 1998a,b) and smoothly extended 
to the low density regime with a gas consisting of neutrons, alpha particles, photons and e ± pairs. 

Under the conditions encountered in the merger neutrinos are emitted copiously and provide the most efficient cooling mechanism for 
the dense, shock- and shear-heated debris. Moreover, the related weak interactions determine the compositional evolution which is altered by 
reactions such as electron and positron captures. The effect on the cooling and the changes in the composition of the material are taken into 
account via a detailed multi-flavour neutrino treatment which is described in detail in Rosswog & Liebendorfer (2003). The results presented 
here are based on the close analysis of late time segments, typically a time i s i m « 15 ms after the start of the merger simulations. Results are 
described from three representative runs: first, an initially corotating system with twice 1.4 M Q (cl.4), second, a system with twice 1.4 M 
and no initial NS spins (il.4), which we regard as the generic case and finally, as an extreme case, a system of twice 2.0 M Q and no initial 
spins (i2.0). In addition, in §4 a set of less resolved simulations is used to explore the dependence of the neutrino luminosity on the system 
mass. 

The debris is heated initially via shocks (for example when two spiral arms merge supersonically, see Fig. 13 in Rosswog & Davies 2002) 
and later via shear motion. The innermost part of the debris torus is heated via shocks (to T « 3 MeV) that emerge when cool, equatorially 
inflowing material collides with matter being shed from the central object. The cool equatorial inflow (T < 1 MeV) allows for the presence 
of heavy nuclei with a mass fraction of ~ 10 % and mass numbers of A ~ 80 and Z/A ~ 0.3. Most parts of the inner disk are essentially 
dissociated into neutrons and protons. We find neutrino luminosities of ~ 10 53 erg s _1 which are dominated by electron anti-neutrinos 
(Rosswog & Liebendorfer 2003). The mean energies are around 8, 15, and 22 MeV for electron neutrinos, electron anti-neutrinos and the 
heavy lepton neutrinos, respectively. Figure 1 shows the total neutrino energy per time and volume in the orbital and meridional planes of 
the merged remnant of our generic case (il.4). The neutrino emission provides the driving stresses that will lead to both a relativistic outflow 
(through annihilation) and a less collimated (and slower) outflow associated with the debris. 



2.1 vv — > e + e above the remnant 

The neutrino annihilation process (which scales roughly with Li) can tap the thermal energy of the hot debris. The presence of a region of 
very low density along the rotation axis, away from the equatorial plane of the debris must be present when the burst takes place or otherwise 
prohibitive baryon loading will happen. The centrifugally evacuated funnel region above the remnant is an attractive place for this deposition 
to occur because it is close to the central energy source (the energy deposition rate scales roughly with the inverse fourth power of the 
distance) but contains only a small number of baryons (e.g. Davies et al. 1994, Ruffert et al. 1997, Rosswog and Davies 2002). 

The ratio of e energy deposition to baryon rest mass energy, r\ — Q V vn n j/{pc 2 ), in the region above the poles of the merged remnant 
is shown in Figure 1. p denotes the matter density and, for simplicity, an energy injection time, ri n j = 1 s has been assumed (Rosswog & 
Ramirez-Ruiz 2002). The copiously emitted electron neutrinos and anti-neutrinos dominate the annihilation process by « 95 %, their mean 
average energies (E) being « 8 MeV for v e and £s 15 MeV for i7 e . 

Due to the finite resolution of the simulations the densities along the rotation axis are overestimated since they are mainly determined 
by particles located in the inner parts of the disk. The above restriction causes r\ to be undervalued. Typical luminosities of ~ 10 47 erg s _1 
(cl.4), fa 3 x 10 48 erg s" 1 (il.4), « 10 49 erg s" 1 (i2.0) are found in the fireball region with r\ > 1. The largest attainable Lorentz factors 
15 although larger values are expected from higher numerical resolution) are found along the binary rotation axis yet at large angles 
away from it an increasing degree of entrainment leads to a drastic decrease in r\ (Fig. 1). The jet is likely to develop a velocity profile so 
that different parts move with different Lorentz factors. During the early stages of its evolution, this entrainment is expected to give rise to a 
mixing instability which could provide erratic changes in the relativistic jet velocity that would manifest in internal shocks. 



2.2 Neutrino-driven wind 

The neutrinos that are emitted from the inner regions of the torus will deposit part of their energy in the outer parts of the thick disk that has 
formed around the central object. For example, a nucleon located at « 100 km from the central object gains the equivalent of its gravitational 
binding energy by capturing 3 neutrinos of « 15 MeV. Similar to the neutrino heating mechanism in core-collapse supernovae, such neutrino 
deposition will drive a strong mass outflow. For high luminosities (L„ > 10 erg s~ ), the mechanical power of the neutrino-driven wind 
L w = 1/2M«J, (where v w is the asymptotic velocity) is found to be largely independent of mass (Thompson et al. 2001) with L w oc 
and ^ < a < Using Thompson et al. (2001) Table 1 we find L w = 2.0 x 10 49 erg s _1 Lj, 6 5 3. This yields a wind mechanical power of 
k 2 x 10 49 erg s~\ « 2 x 10 50 erg s~\ and w 2 x 10 51 erg s _1 for runs cl.4, il.4 and i2.0. 

While the above estimates agree with other estimates derived in the context of proto-neutron star formation (e.g. Duncan, Shapiro & 
Wassermann 1986; Qian & Woosley 1996) they do not take into account the specific geometry of a merger remnant nor the extreme neutron 
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Figure 1. Geometry of the neutrino emission above the merged remnant of model il.4 (2x 1.4 Mq, no spin). The right panel shows the logarithm of the total 
neutrino energy emitted per time per volume (cgs units) from the orbital plane while the left panel displays the meridional plane. The upper half-plane of the 
left panel shows the ratio of energy deposited during T; n j via vv — > e+e - to rest mass energy, rj, which is a measure of the maximum attainable Lorentz 
factor. The colourbar on the left of this panel refers to log(»y). The lower half-plane of the left panel shows the total neutrino emission in meridional plane. The 
colour coding is the same as in the right panel. 

richness of the debris, we expect the corresponding correction factors to be of order unity, which we consider accurate enough for the 
following discussion. 



3 THE JET- WIND CONNECTION 

We have now all the ingredients to consider the interaction of the jet with the less collimated and slower outflow associated with the debris. 
The most prominent feature of this interaction is that the wind may be responsible for collimating the jet (LE). The collision of the two fluids 
will lead to the formation of a contact discontinuity across which the total pressure is continuous. Besides, two oblique shocks, one in each 
fluid, will form across where the streamlines of the colliding (unshocked) fluids are deflected. The structure of the shocked layers will depend 
on the parameters of the two outflows and on the boundary conditions. 

On the condition that the baryonic wind is highly supersonic not too far out the streamlines of the colliding fluids, the momentum 
transfer by the wind into the jet is dominated by ram pressure. The ram pressure of the baryonic outflow just upstream the oblique shock at 
the axial coordinate z is then given by (LE) 

(pw + p-wc)Vt = f3 w L w /(A-K 2 ca c r), (1) 

where p w , p w c 2 , and U w the corresponding wind pressure, energy and four- velocity, r is the distance between the torus and the nearest point 
at the shock, a c (z) is the cross-sectional radius of the contact discontinuity surface, and f3 w is the wind terminal velocity (which is set to the 
escape velocity of the system in the following estimates). 

If the shocked wind layer is very thin due to rapid cooling of the gas and using the requirement that the total pressure must be continuous 
across the contact discontinuity, LE balance entropy and energy fluxes to obtain: 9 ~ 7T/3" 1 (Lj / L w ), where 8 is the jet opening angle at 
large distances and the Lorentz factor at the base of the jet has been assumed to be of order unity (this approximation is justified in Rosswog 
& Ramirez-Ruiz 2002). 

In the above estimate L w is the total wind power emanating from a thin torus of radius R centred around the central object. This 
geometry is reminiscent of the neutrino emission structure seen in Figure 1 . Neutrinos emitted from the most luminous ring-like region at 
about R ~ 25 km drive a baryonic outflow whose total wind power surpass that ejected from rest of the debris by at least an order of 
magnitude. Numerical integration of equation Q, where the neutrino-driven wind is assumed to arise from a sequence of concentric rings 
located at z = 0, confirms that the ring with the highest wind power determines the properties of the cross-sectional radius of the jet at large 
distances (z > R). To first order we therefore estimate L w to be the total wind power emanating from this thin and most luminous region. In 
this way, at large enough distances from the plane of the remnant, the semi-aperture angle of the wind-confined jet is given by 

9 W 0.3 Z^, 1 -! L W * B0 ij,48.5, (2) 

where we adopt the convention Q = 10 x Q x , using cgs units. Once the fireball Lorentz factor exceeds 9~ 1 , it will remain conical with the 
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Figure 2. Total neutrino luminosities as functions of time for various NS/NS merger masses. The merger calculations were performed using a SPH method 
with ~ 6 X 10 4 particles. 

same opening angle regardless of the external conditions. This occurs at a distance ~ RQ~ X from the engine, which in turn implies that the 
baryonic outflow needs to develop t ~ 2>R\f5Z\9Z\ ms before the explosion that forms the GRB. This requirement is justified by the fact 
that the neutrino luminosities reached their maximum, stationary level t > 15 ms after the start of the simulations (see Fig. 2). 

It is therefore likely that the neutrino-driven outflow that develops in the merged debris will have enough pressure or inertia to provide 
collimation. The beaming fraction is given by ? = Q/2tt ~ (iv 2 /2)(Lj / f3 w L w ) 2 and hence the apparent luminosity is Ln ~ Lj^^ 1 = 
2 x 10 49 /3 2 _il£,5o Ljl s erg s" 1 . The above estimate gives « 0.1 (» 0.02) and Ln ~ 10 51 erg s" 1 (Ln ~ 3 x 10 52 erg s" 1 ) for runs 
il.4 and i2.0 (assuming a — 16/5), which clearly satisfies the apparent isotropized energies of w 10 51 erg implied for short bursts at z = 1 
(Panaitescu et al. 2001; Lazzati et al. 2001). 

The above reasoning does not take into account effects associated with the formation of a rarefraction wave in the baryonic wind near 
the interface separating the two fluids. This will depend on the properties of the solution near the object and could result in a steeper decline 
of the pressure supporting the jet, and the consequent alteration of the profile at the contact discontinuity. There will also surely be some 
entrainment of the electron-ion fireball plasma and this should ultimately show up in the polarization observations which can, in principle 
distinguish a pair plasma from a protonic plasma. Entrainment will be promoted by linear instabilities that could grow at the jet surface 
(Appl & Camenzind 1993). In addition, there can be an exchange of linear momentum with the surrounding outflow, even if there is minimal 
exchange of mass. This can develop a velocity profile in the jet so that different parts move with different Lorentz factors. 



4 DIVERSITY OF AFTERGLOW BEHAVIOUR 

The framework we have used to determine the geometry and energetics of short GRBs is a simple one. Given the limitations of our assump- 
tions we draw three conclusions. First, an inner relativistic jet collimated by a surrounding baryonic wind can appear brighter to the observer 
by a factor that is inversely proportional to its intrinsic power (LE). Second, a modestly broad distribution of neutrino luminosities can result 
in a wide variety of opening angles and apparent luminosities. This can be understood as follows: the luminosity of the ^-annihilation jet 
increases as L 2 , while L w oc L™, so that 6 oc L"~ 2 and Ln oc L 2a ~ 2 (^ < a < ¥-). Third, only a small fraction of bursts are visible to a 
given observer with most of the energy of the event being in the unseen outer outflow. 

Figure 2 shows the results of a set of test calculations (with ~ 6 x 10 4 particles, no initial spin) performed in order to explore the 
effect that the masses of the neutron stars in the merging binary have on the total neutrino luminosity. From these models we derived an 
approximate relation between L„ and the total mass of the compact binary, which we then apply to estimate the possible variety of afterglow 
behaviour that could be expected from NS-NS mergers. 

We set up a 10 6 grid of initial binary masses chosen according to the theoretical distribution function of NS derived by Fryer & 
Kalogera (2001). The compact-remnant distribution is dominated by NSs in the mass range 1.2-1.6M and falls off drastically at higher 
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Figure 3. Isotropic luminosities implied for short GRBs understood in terms of collimation of the fireball by an outer baryonic outfiow.The left panel shows the 
distribution of neutron-star remnant masses derived by Fryer & Kalogera (2001). The NS masses are measured by Thorsett & Chakrabarty (1999) and Barzfv 

et al. (2001). The right panel shows the distribution of expected opening angles and luminosities for NS binaries. The dashed line assumes L w oc while 

17/5 

the dotted line uses L w oc L v . 



remnant masses (Fig. 3). Under the assumption that the total neutrino luminosity (at late times) depends primarily on the total mass of the 
compact binary (note that this is a reasonable assumption for mergers of unequal masses provided that the binary mass ratio is not extreme; 
see Fig. 2) and that L w oc L" (Thompson et al. 2001), we derive the distribution of opening angles and luminosities expected from such 
encounters. Our results for two values of a (16/5 and 17/5) are shown in Figure 3. 

Since the jets are visible to only a fraction <, of the observers, the true GRB rate R T — {<;}~ 1 Robs, where i? b s is the observed rate and 
is the harmonic mean of the beaming fractions. We find ~ 100. Following Schmidt (2001) i? bs(z = 0) = 0.5 Gpc~ 3 yr -1 . 

The true rate is then R T (z = 0) ~ 50 Gpc~ 3 yr -1 , which should be compared with the estimated rate of NS coalescence i?Ns(z = 0) ~ 80 
Gpc -3 yr -1 (~ 2 x 10~ 5 yr -1 per galaxy; Phinney 1991). We note that the simple collimation mechanism presented here provides an 
adequate description of the observations. 



5 SUMMARY 

In this letter we consider the viability of NS-NS binary coalescences as central engines of short GRBs. The main form of energy release we 
regard is neutrino emission from the disk 1 . 

We suggest that the isotropic luminosities implied for short GRBs at cosmological distances can be understood in terms of collimation 
of the ^-annihilation fireball by the neutrino-driven wind. We argue that the existence of such a baryonic wind is a natural consequence of 
the huge gravitational binding energy released in the form of neutrinos. 

Moreover, we show that the wind derived from the disk could in principle have enough pressure or inertia to provide collimation to 
the fireball that derives from the region around the BH. In other words, bursts produced by such a mechanism are not isotropic but instead 
beamed within a solid angle, typically ~ 0.1 sterad. Within the framework of this simple model we have deduced the distribution of opening 
angles and luminosities expected from ^-annihilation in NS-NS mergers and uncovered that within the uncertainties of the calculations, the 
coalescence scenario is capable of providing the required isotropized energies and a compatible number of progenitors. 

Given the few hundred kilometres per second acquired by the NS at birth and assuming a typical coalescence time of about 100 Myr, 
these bursts should occur predominantly in the low-density halo of the galaxy. The external shock can therefore occur at much larger radii 
and over a much longer timescale than in usual afterglows, and the X-ray intensity is below the threshold for triggering. Alternatively, the 
burst could go off inside a pulsar cavity inflated by one of the neutron stars in the precursor binary. Such cavities can be as large as fractions 
of a parsec or more, giving rise to a deceleration shock months after the GRB with a consequently much lower brightness that could avoid 
triggering and detection. In conclusion, the absence of detected afterglows in short bursts is not surprising, and, as argued above, a wide 
diversity of behaviours may be the rule, rather than the exception. 



1 More energy could be pumped into the e ± 7 fireball when the rapidly rotating BH is formed or if magnetics field are able to tap the rotational energy of the 
BH with higher efficiency than vv does (Popham, Woosley & Fryer 1999; Rosswog, Ramirez-Ruiz & Davies 2003). 
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